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Summary. It is problematic to interpret the quantum jumps of an atom inter- 
acting with thermal light in terms of counts at detectors monitoring the atom's 
inputs and outputs. As an alternative, we develop an interpretation based on a 
self-consistency argument. We include one mode of the thermal field in the system 
Hamiltonian and describe its interaction with the atom by an entangled quantum 
state while assuming that the other modes induce quantum jumps in the usual 
fashion. In the weak-coupling limit, the photon number expectation of the selected 
mode is also seen to execute quantum jumps, although more generally, for stronger 
coupling, Rabi oscillations are observed; the equilibrium photon number distribu- 
tion is a Bose-Einstein distribution. Each mode may be viewed in isolation in a 
similar fashion, and summing over their weak-coupling jump rates returns the net 
jump rates for the atom assumed at the outset. 



1 Introduction 

The notion of a quantum jump entered physics with Bohr's model of the atom and 
was elaborated in a semi-quantitative formulation in Einstein A and B theory JlJ . It 
was, from the beginning, an idea at variance with the usual commitment to a con- 
tinuous time evolution and with the continuous constitution of light as an electro- 
magnetic wave. The appearance of the Schrodinger equation relieved the situation 
somewhat, but ultimately, through the use of perturbation theory to make testable 
predictions about quantum scattering processes, the quantum jump remains with 
us, though certainly in a more sophisticated and mathematically refined form. 

The Monte-Carlo wavefunction and quantum trajectory methods developed in 
quantum optics use quantum jumps as an explicit component of a stochastic 

time evolution in a manner very reminiscent of Einstein A and B theory. In fact, the 
only novelty is to combine Einstein's rules for quantum jumps with a coherent evo- 
lution between jumps that admits a dynamic involving superpositions of stationary 
states. In this, these methods achieve something remarkably similar to the proposal 
of Bohr, Kramers, and Slater (BKS) |,| for uniting discontinuous jumps among 
the stationary states of a material system with a continuous evolution between 
jumps, during which time material oscillators, possessing coherent amplitudes, are 
brought into play. 

The realistic interpretation sought by BKS may not, however, be entertained. 
In most circumstances, an interpretation of the jumps employed in quantum trajec- 
tory theory is based upon a record of time-resolved photon counts which might be 
realized in practice by terminating every output channel in a photodetector PJl5|] . 
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This scenario is plausible because optical frequencies are sufficiently high that scat- 
tered photons can be detected against an essentially vacuum-state background. The 
measurement-based interpretation is problematic, though, for an atom exchanging 
photons with a thermal environment. In this situation, incoherent photons are both 
emitted and absorbed; moreover, it is impossible to distinguish a scattered photon 
from some other photon in the environment. Of course, schemes such as electron 
shelving exist that are able to monitor thermal quantum jumps pp|jic|] . They, how- 
ever, make intrusive measurements by utilizing strong couplings to other inputs and 
outputs, and are not a suitable foundation for the interpretation of quantum tra- 
jectory equations. The relationship, in fact, is exactly the reverse; electron shelving 
is one of the measurement schemes that quantum trajectories would propose to 
explain. 

In this paper we follow a different direction to give substance, beyond a mere 
assertion, to the interpretation that an atom exchanging photons with a thermal 
environment does, in some well-defined sense, execute jumps. We define the sense 
through a self-consistency argument. Contrary to the notion of quantum jumps, 
quantum mechanics continuously entangles two interacting systems through the 
Schrodinger evolution. We show that a single mode, selected from the many modes 
of a thermal environment, and allowed to evolve in interacting with an atom to 
produce such an entanglement, is, in fact, seen to undergo jumps in its photon 
number expectation if it is assumed that all other modes of the environment, treated 
collectively as a reservoir, induce jumps between the atomic states according to the 
rules of Einstein theory. The jump evolution for the selected mode emerges from an 
otherwise continuous evolution in the weak-coupling limit. Thus, assuming jumps 
induced by the reservoir as a whole leads, self-consistently, to the appearance of 
jumps in an individual mode of the reservoir when that mode is allowed to entangle 
with the atom via the standard interaction Hamiltonian and Schrodinger evolution. 
The jumps do bring the individual mode to a Bose- Einstein distribution over photon 
number, and the jump rates, derived for every mode viewed individually in this way, 
sum to net rates which agree with Einstein A and B theory. 

The underlying theme of this paper is the conflict between a continuous and 
a discontinuous quantum evolution and the self-consistency of the two in the per- 
turbative weak-coupling limit. We therefore briefly review, in Sect. ^, the quan- 
tum jump model of Einstein, the BKS proposal to include a continuous evolution, 
and the quantum trajectory realization of the latter. The argument for the self- 
consistency of thermal jump trajectories is elaborated in Sect. ra. 



2 Einstein A and B Theory, BKS, and Quantum 
Trajectories 

2.1 Thermal Quantum Jumps 

We consider the single two-state atom illustrated in Fig. [l], in thermal equilibrium 
with Planck radiation at temperature T. In Einstein A and B theory photons are 
exchanged between the atom and the radiation field as the atom jumps randomly 
between its two stationary states. The jump rates follow a prescription taking into 
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/down = A + Ba(uj Q ) 
T up = Ba (w ) , 



(la) 
(lb) 



where 



r(u>o) = n{uj )huj Q [p{u]o) /V] 



(2) 



is the energy density of the radiation field at the resonance frequency ojo of the 
atom, with average photon number per mode 



n(w ) = [e^o/kBT _ -^-l 



(3) 



and mode density (in volume V) 



(4) 



and the Einstein A and B coefficients must satisfy 



B TT 2 C 3 



A Tiul 



(5) 



in order for the atom to be brought into thermal equilibrium with the radiation. 




Fig. 1. Thermal quantum jumps in Einstein A and B 
theory. The jump rates are defined in (^) and ([lb]). 



With the help of the relationship (|5|), we write ( |la| ) and ( |lb| ) in the modern 
notation 



-Tdown = A[n(uo) + 1] 

r up = An(u)o) . 



(6a) 
(6b) 
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Einstein theory does not assign a value to the coefficient A. From quantum me- 
chanics, however, using Fermi's Golden rule we obtain 



(7) 



where 



|«A,*l(wo)| 



2ite V lex - A esl 



(8) 



is the dipole coupling strength to a mode of the radiation field with polarization 
A and direction of propagation specified by the unit vector n (polarization vector 



d eg is the atomic dipole matrix element. 



Commonly, Einstein theory is discussed at the level of rate equations for the 
occupation probabilities of the atomic stationary states. The theory does, however, 
define a stochastic process - one that may be visualized in terms of quantum jumps 
whose occurrences unfold randomly in time. With each realization of the stochastic 
process we associate a record of jump types and jump times, 



REC 



r d , 

t + Tl 



t + T 2 t + T 3 



up -'down 

t + t 4 t + r s 



(9) 



In Fig. |^ we illustrate the discontinuous evolution of the atom in coordination with 
its absorption and emission of thermal photons. 



r, 



r 



r 



t+T, t+T, 



Fig. 2. Sample realization 
of the Einstein stochastic 
process. The corresponding 
record is defined in (pi). 



2.2 Coherence: The BKS Proposal 

Although Bohr had himself put forward the notion of a quantum jump to explain 
the association of stationary-state energy differences with electromagnetic wave 
frequencies in his model of the hydrogen atom, he became quite dissatisfied with 
the idea in the concrete form it acquired under Einstein's proposal. It was the 
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light quantum, specifically, that troubled him the most. Bohr insisted that since 
so many optical phenomena rely on the continuity of coherent waves, the wave 
nature of light simply could not be dismissed. He recognized, on the other hand, 
that a discontinuous process was definitely needed to account for light emission 
and detection. What he was unavoidably drawn towards, then, was some sort of 
merging of the two ideas. 

A program to accomplish this was outlined in what has come to known as 
the BKS proposal J5|,[|. The central idea of this proposal is the proposition that 
during the residence times in a stationary state, represented by the horizontal lines 
in Fig. ^, an atom is not inactive in its interaction with the electromagnetic field; 
rather, it acts through a coherent dipole radiator, or "virtual oscillator" in the words 
of BKS, which is all the time radiating an electromagnetic wave of frequency uiq = 
(E e — E g )/h. This wave is either in phase or out of phase with the external radiation 
at frequency ujq depending on whether the residence is in the stationary state with 
energy E e or E g . Thus, there is an energy transfer under the laws of classical 
electrodynamics either to the electromagnetic field from the dipole, or in the reverse 
direction, depending on the stationary state JllJ . Clearly a double counting of the 
exchanged energy occurs if one accepts that light quanta are also emitted and 
absorbed at the times of the jumps. However, the goal was precisely to eliminate 
these quanta, although still permitting the atom to jump. BKS attempted, thus, 
to retain, but keep separate, two incompatible mechanisms for energy exchange - 
a wave (continuous) mechanism for the absorption and emission of radiation and a 
particle (discontinuous) mechanism for the change of material state energies. Their 
proposal foundered on its obvious violation of energy conservation at the level of 
the individual quantum events, a feature that appeared not to be supported in 
Compton scattering experiments E,5j[. 



2.3 Coherence: Quantum Trajectories 

In retrospect we can see that Bohr had in mind a conception of light that, although 
supported by numerous wave phenomena in optics, was largely inappropriate for 
the light sources available at the time. Radiators in thermal equilibrium are not 
sources of coherent waves. They radiate electromagnetic noise, which, with filtering, 
can approximate low intensity light possessing first-order coherence, but is very far 
from the classical concept of a coherent wave of large and adjustable amplitude. 
Modern lasers, however, emit something close to the classical ideal. In their case, 
the high coherence, if it is to be preserved, disallows the tracking of energy at the 
level of the individual quanta, so that the energy conservation argument against 
the BKS proposal does not apply. 

Quantum trajectory theory is designed to deal with problems involving the 
interaction of matter with the high coherence light sources available in modern 
laboratories. It shows remarkable similarities to the BKS proposal. These are de- 
scribed elsewhere , and we do not plan to discuss them in any depth here. As 
an introduction to our main topic, however, it is useful to contrast Einstein A and 
B theory with the quantum trajectory description of a coherent field, amplitude £, 
resonantly exciting the two-state atom of Fig. |l| The connections with BKS emerge 
automatically through this exercise. 

The atom is still located in a thermal environment and quantum jumps still 
appear as they do in Fig. 0- but with one notable modification. Due to the induced 
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coherence, it is necessary that the system state be a superposition of the stationary 
states \E e ) and \E g ), which we denote by [V'RBc(t))- As suggested by BKS, there is 
a coherent interaction with the electromagnetic field between the quantum jumps. 
This we account for by a continuous evolution under the Schrodinger equation 
|s]|3|.[l| (for the unnormalized conditional state) 

T7 = -T-Hb \VREC) , (10) 

at in 

with non-Hermitian Hamiltonian 

H B = \%{uo - ir dawn )\Ee)(E B \ - \h(uo+ir up )\Eg){Eg\ 

+th£(e^ ot \E g ){E e \ - e- luJat \E e )(E g \) , (11) 



in which the external coherent field is classical, and its interaction with the atom is 
treated in the dipole and rotating-wave approximations. The quantum jumps are 
governed by the probabilistic rules of Einstein A and B theory, generalized, in a 
natural way, to account for the fact that the system at any time is not definitely in 
a particular stationary state. There are jumps 



I^rec) ^ (\E g )(E e \)\j wc ) , (12a) 
IV'rec) ^ (\E e )(E g \)\i> REC ) , (12b) 



with jump rates 

-Rdown = -Tdown|(-Ee|V'REc)| 2 , (13a) 
Rup = r u p|(E g |V>REC>| 2 . (13b) 

Figure ^ illustrates how realizations of this stochastic process appear. In (a) 
the coherent excitation is relatively weak and the overall form of the evolution re- 
mains close to that produced by the Bohr-Einstein quantum jumps (Fig. ^). There 
is, however, in addition to the switching of the energy, a weak induced coherence 
carried along by the continuous evolution between jumps as a nonvanishing polar- 
ization amplitude, quite reminiscent of the BKS virtual oscillator. In Fig. Ea(b) the 
coherent excitation is much stronger. Here, the dominant mechanism for evolution 
between the stationary states is a coherent Rabi oscillation; we move to a coher- 
ent evolution that is nonperturbative and a regime which only became accessible 
with the invention of the laser. This, specifically, is a form of evolution following 
from the Schrodinger equation; nonperturbative coherence was not anticipated by 
the BKS proposal. We might note in passing that lasing without inversion and re- 
lated phenomena acquire their counterintuitive features from such nonperturbative 
coherence 
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(a) (b) 




Fig. 3. Sample quantum trajectories with both thermal jumps and induced coher- 
ence: n(w ) = 0.25, £/A = 0.1 and 1.5 [(a) and (b)]. 



3 Self-Consistency of Thermal Quantum Jumps 

At a sufficiently low temperature, when n(u;o) -C 1, the jump record that labels the 
state can reasonably be made by detectors monitoring the scattered light. Almost 
all jumps will be down-jumps governed by the spontaneous emission rate A which 
may be identified with emitted photons counted as isolated (in space and time) 
excitations of the vacuum. For many applications in quantum optics this is the 
situation in reality. Nonetheless, the thermal jumps, though they might be negligible 
in practice, cannot be set aside from a fundamental point of view. Considering then 
those jumps which cannot reasonably be identified with the "click" of a detector, 
is there any other justification for returning to the language of the old quantum 
theory, given that, in quantum mechanics, the Schrodinger equation invokes only 
a continuous evolution? We aim to show that there is, in so far as the jumps are 
self-consistent - consistent with the Schrodinger evolution - in the weak-coupling 
limit. 

3.1 Trajectories for a Single Field Mode 

The Hamiltonian for a two-state atom interacting with the radiation field of a 
thermal environment (we now take £ = 0) is 

H = ±hbj \E e )(E e \ - \Ku \E g ){E g \ 

A' ,n' 
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+ R[«A',A'(w')|S e >(S 1 ,|f v , A / i(U ,+H.c.], (14) 



where f\ - and f^^,^ are creation and annihilation operators for the field mode 
with polarization A, propagation direction n, and frequency u, and Kx,n(w) is the 
mode coup ling coefficient whose magnitude is defined in (|^). The stochastic process 
(|ic[)— ( |l3b ) (£ = 0) is developed, formally, around the master equation derived 
from (|14D • This master equation describes the quantum state of the atom 

alone, after tracing over every mode of the radiation field. Our idea is to raise one 
mode of the field to the same status as the atom by including it, along with its 
interaction with the atom, in the system Hamiltonian. All other modes are to be 
treated as a reservoir as before and their interaction with the atom described by 
quantum jumps. The stochastic process is the same as in (]ic|)-( p^b| ), but with the 
non-Hermitian Hamiltonian Hb replaced by 

H B = \%(u - ir dawn )\E e )(E e \ - ±h(u +ir np )\E g )(E g \ 

+h[K X , A (u;)\E e ){E g \fx,&, u +H.c.] . (15) 

Of course removing one mode from the reservoir has no effect on the overall jump 
rates for the atom. The change is that we can now follow the evolution of an explicit 
Hilbert space vector for the selected mode, one that entangles this mode with the 
atom. We ask how does the selected mode evolve in the Hilbert space; in particular, 
does it also experience quantum jumps? 



Figures 



and |B| show sample trajectories for the selected-mode photon number 
expectation [16| - for a series of decreasing coupling strengths, (a)-(d), and assum- 
ing resonance with the atom, Fig. ^, and a detuning from the atom, Fig. ^| With the 
coupling strong compared to the Einstein A coefficient coherent Rabi oscillations 
are seen. There are also discontinuous changes, which in the case of strong coupling 
are merely a direct manifestation of the assumed quantum jumps for the atom. 
Note, however, that the atom does not jump monotonously, "up" then "down" 
then "up" ■ • ■ , as in Fig. ^; repeated up-jumps can transfer many energy quanta 
to the field mode. At an intermediate coupling strength, partial Rabi oscillations 
are still present. Once the coupling becomes weak, though, the Rabi oscillations 
apparently disappear altogether, and an entirely new kind of jump evolution sets 
in. These jumps proceed at a rate far less than the overall jump rate for the atom. 
Their rate decreases with the square of the coupling constant [(c) to (d)] and also 
when the detuning is increased (from Fig. ^ to Fig. |^). 

Having, then, assumed jumps for the atom in interaction with all but one of 
the field modes, a jump evolution for the one remaining mode emerges naturally 
in the weak-coupling limit. We close the self-consistent loop by showing that the 
one mode samples a Bose-Einstein distribution, and by calculating the rate of the 
single-mode jumps, to demonstrate that the sum over modes returns the rates T^p 
and r do wn. 
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Fig. 4. Sample trajectories for the 
photon number expectation of a 
single field mode included as part 
of the system. The mode is reso- 
nant: Aui/A — 0, n(ojo) = 1, and 
\k x , A (uo)\/A = 10, 1.0, 0.1, and 
0.01 [(a), (b), (c), and (d)]. 



3.2 Self-Consistency 

Let us denote the number of energy quanta shared between the atom and the field 
mode at time t by n„ + 1. Let tk be the time of the very last jump of the atom and 
tk+i be the time of the jump that is to occur next. Then, for tk < t < tk+i, the 
entangled state of the atom and field mode may be expanded as 



I^rbcW) = C eli {t)\E e )\n u ) + C g \i(t)\E g )\n u + 1) 



(16) 
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Fig. 5. Sample trajectories for the 
photon number expectation of a 
single field mode included as part 
of the system. The mode is non- 
resonant: Aui/A = 2.0, n(uo) = 1, 
and \K X ,ri{u)\/A = 10, 1.0, 0.1, 
and 0.01 [(a), (b), (c), and (d)]. 



with C e |i(0) = <5 e ,i and C 9 |i(0) = 5 g< i, where i is e or g for an up- or down-jump at 
tfc, respectively. Writing 



(17a) 
(17b) 
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with 4>a{w) = ar g[ K A,ft( w )]; from ( |lo| ) and (p^), the equations of motion for the 
conditional state amplitudes are 



—rr- = -|(Aown - iAjJCeii - i\it\ t A (o>) \VriZ + 1 , (18a) 



dCgU 



dt 



f 11 = -|(r up + iZia;)C , 9 | l -ii KA , A (a;)|V^TTC' e | i . (18b) 



Consider now the case of an up-jump at tfc, such that the initial state amplitudes 
are C e u(Q) = 1 and C g \i(0) = 0. For weak coupling we will have C e u(t) w 1 and 
^gli(^) ~ |ka,a(w)|, with an overwhelming probability that the next jump of the 
atom will be a down-jump. By a similar argument, it is highly likely that the down- 
jump is followed by another up-jump; the most probable progression is then "up", 
"down", "up", . . ., just as we illustrated it in Fig. (notice that the number is 
unchanged throughout such a progression). Due, however, to the small amplitude 
- either C g \i(t) ~ |«a,si(w)| or C e \i(t) ~ |ka,a(w)| - excited by the coupling of 
the atom to the selected mode, there is always a small probability that a jump 
will occur to break the alternating sequence. An up-jump might be followed by a 
second up-jump or two down-jumps might occur in a row. These events change 
and produce the jumps of the field mode seen in Figs. ^ and ^. In Figs, ^(c) and 
^(c), the presence of the small amplitude that underlies the jump mechanism is still 
seen as a "fuzz" on top of the developing smooth curve. The "fuzz" is even present 
in Figs, ^d) and ||(d), though there it is too small to be visible. 

The physical interpretation of the anomalous events in the jump record of the 
atom is that each represents the scattering of a photon between one of the many 
field modes of the reservoir and the field mode selected to be viewed. Two up- 
jumps occur in a row for example, because, in the interval between them, the 
energy absorbed on the first jump is transferred to the selected mode; the quantum 
trajectory resolves the transfer at the time of the second up-jump. 

Our task now is to calculate rates for the unlikely j ump s. We do t his using the 



method of Sect. IVD in p.4fl . The equations of motion (18a) and (18b) give 
d\C e \i\ 2 



dt 
d\C g]i \ 2 
dt 

dRe(C eli C* gli ) 
dt 

dt 



-r down \C ell \ 2 - 2| KA , A (o;)|\/^TTlm(C , eU C; u ) , (19a) 



-r up \C gl ,\ 2 + 2| KA , A (^)| % AUTTlm(C e | l C*;| l ) , (19b) 

-^(r down + r np )Re(C el ,C* lt ) - AtuIm{C e]l C^) , (19c) 
-^(r down + r^ p )lm(C e{i C g]i ) + ALjRe{C e]l C* g}l ) 

+K*iMlx/^TT(|c e |,| 2 - \c g]t \ 2 ) . (i9d) 
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We define 



W e = I dt\C ell (t)\ , W a = / dt\C 9li (t)\\ 



(20) 



and 



U ee Re 



dt C,, \i(t) Cg I i 



V = Im 



dtC e |j(t)C s |j 



(21) 



where -TupWg and /down We are the probabilities, given i is e and g , res pect ively , 
that the unlikely jump will occur. From the Laplace transforms of (19a) - (19d), 
we then have 



■S e ,i = -r dovjn W e - 2\n x , fl {uj)Wn UJ + IV , (22a) 

-<5 9ll = -r up W g +2\K X , A (fij)\y/riZ + 1 V , (22b) 

o = -i(r down + r up )v + AluU + I^HIV^ +T(w e - w g ) , (22c) 

o = -±(r down + r up )[/-zw, (22d) 



and hence 



^upW s | j=e = 
i~downWe ,•_„ = 



2 (/down + r up )/ir 2TY\K\,fi(uj)\ 2 
[§(/daw„ + r»p)] 2 + (A^) 2 ^own 

[|(r d0 wn + r„ p )] 2 + (z\o;) 2 r " P 



(n w + 1) , 
(n w + 1) , 



(23a) 
(23b) 



where we have solv ed (|22a|)-(22d ) to lowest order in the coupling strength. 



Equations (23a) and (|23b|) specify the probability for the unlikely jump to occur 



following any preparation of the initial state i. To obtain photon-number jump rates, 
we must multiply b y th e rate at which the state i is pre pared, i.e., by the jump 
rates for the atom; (23e) is multiplied by _r up p° q and (23b) by -TdownP^ 1 , where p g q 



and Pe q are t he st ate occupation probabili ties i n thermal equilibrium. We also set 
n w — N u in (23a) and n u + 1 = N u in (23b), where iV w is the photon number 
expectation plotted in Figs. ^ and 1 (recalTtnat n u + 1 is the number of quanta 
shared with the atom at tk)- The photon-number jump rates are then 



up 2 (-^down + r up )/lT \|2/»t . -,\ cq / n . \ 

7a£ = — " - 2 2it\k x , a (uj)\ (N u + l)p e q , (24a) 

[i(Aown+r up )] + (Au) 2 

down 2 (^ dow ° ^ u p)/ n r, 1 I m 2 Ar cq /„., > 

7iV„ = ri , l2 27r|K AiTi (a;)| Ar„p 9 4 . (24b) 

\ 2 



[|(r d own + r up )] + (ZVu/p 
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The self-consisten cy of the t herm al jump picture is now easy to demonstrate. On 
the one hand, we use ( |24a| ) and ( |24b| ) to set up rate equations for the selected mode 
photon number and, using detailed balance, solve these in steady state. Hence, we 
obtain the equilibrium probability to find N u photons in the selected mode: 



= (l'PtVp?)[%) =[nM + l] 



n(wo) 



n(uj ) + 1 



(25) 



where we have used pt q /pg q 



/., — -Tup/rdown and the Einstein formulas ( |6aj ) and (jsbj). 
We obtain a Bose-Einstein distribution with average photon number 7V„ = n(u>o); 
note that the specific mode coupling strength and frequency affects only the rate 
of approach to equilibrium. Of course, in reality, each field mode couples to a 
vast number of two-state systems, and most strongly to those with which it is 
nearly resonant. For the realistic situation we would therefore find the expected 
Nu = n(u), consistent with the Planck radiation formula. 

On the other side we must dem onstrate the self-consistency of the jump rates. 
To this end, we sum (24a) and (24b) over all modes (all A, n, u), with N u replaced 
by its average value, and neglecting the frequency dependence of the density of 
states and dipole coupling constant (in light of the Lorentzian resonance). The 
resulting jump rates for the gain and loss of photons by the thermal environment 
should equal the jump rates assumed initially for the atom. The sums do, indeed, 
return / down p cq and F up p^, showing that the net jump rates are in accord with the 
Einstein rules (|&j) and ( jog ) and Fermi's golden rule, ((?]) and (|^). This completes 
our demonstration that thermal jump trajectories are consistent. 
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